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ABSTRACT 

Wc investigate the capability of pulsar timing arrays (PTAs) as a probe of primordial 
black holes (PBHs), which might constitute the Galactic dark matter. A PBH passing 
nearby the Earth or a pulsar gives an impulse acceleration and induces residuals on 
otherwise orderly pulsar timing data. We show that the timing residuals induced at 
pulsars are optimal for searching heavier PBHs than those at the Earth, and the two 
probes are highly complemcntal. Future facilities like SKA could detect PBHs with 
masses around ~ 10 22_28 g even if only a small fraction (< 1%) of the Galactic dark 
matter consists of these PBHs. 
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1 INTRODUCTION 

A significant fraction of matter in our Galaxy is consid- 
ered to be occupied by dark matter, but its nature is p oorly 
understood at present (|Bertone. Hooper, fc Silk 2005). Pri- 
mordial black holes (PBHs) are an interesting astrophys- 
ical candidate of dark matter, and various observational 
constraints have been posed on their allowed mass range 
Oarr et alj|20ld : iKhlopovllioiol ) . However, the current con- 
straints in the mass range 10 20 g < Mpbh < 10 27 g remain 
relatively weak. 

To examine PBHs in this range, ISeto fc Cooravl (|2007T ) 
proposed an observational technique to directly probe grav- 
itational interaction between the solar system and a nearby 
PBH. Around their close approach, the Earth receives an 
impulse of acceleration whose profile depends on the mass, 
distance, and velocity of the PBH. Given the local mass 
density of dark matter, the expected magnitude of the ac- 
celeration is very small, but high-precision measurements 
of pulsar timing could allow us t o det ec t the weak impulse 
signal (see also Seto fc Cooravl |2004 ) ; ISaito fc Yokovamal 
l|2009D : lGrTest et all |201ll )). 

Roughly speaking, in the present context, the pulsar 
timing analysis can be essentially regarded as measurement 
of the arrival times of radio pulses that were emitted by a 
pulsar and received on the Earth. The local accelerations of 
both of the pulsar and the Earth are encoded in the modu- 
lation of the time-of-arrival (TOA) data, as a simple linear 
combination of two separate terms for the two masses. We 
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call them by the pulsar term and the Earth term respec- 
tively. 

When observing multiple pulsars for PBH search, the 
Earth terms are commonly excited by the acceleration of 
the Earth, and have coherent structure among TOA data of 
different pulsars. Therefore, we can statistically amplify the 
weak impulse signature on the Earth by using a pulsar tim- 
ing array (PTA) and effectively reducing the timing noises. 
The underlying statistical approach here is similar to that 
for detecting gravitational waves (GWs) whose effect on the 
TOA data can be also expressed with two terms induced sep- 
arately at the pulsar and the Earth <|Sazhinlll978l ; iDetweileil 
ll979l : (Hellings fc Downs!ll983D . For the GW signals, we call 
them by the pulsar GW term and the Earth GW term, in 
order to distinguish them from the two acceleration terms. 
The traditional target of PTAs is the Earth GW terms. 

Recently, it has been actively discussed that we 
might utilize the pulsar GW terms for analyzing GW 
sources such as merging super-massive black hole binaries 
Jjenet et alj|2004 ICorbin fc Cornishll2010l ; iLee et alJlioTll : 
lEllis. Jenet. fc McLaughlirj|2012l )7 To this end, we need suf- 
ficiently stable millisecond pulsars (MSPs), and such pre- 
ferred systems might be discovered with future observational 
facilities. With the pulsar GW terms and known distances 
to the individual pulsars, our information content on GWs 
at the nano-Hertz band can be greatly increased, compared 
with that extracted from the Earth GW terms alone. 

Meanwhile, if the timing noises of individual pulsars are 
small, we can also probe PBHs around the pulsars through 
the pulsar (acceleration) terms, in addition to PBHs close to 
the solar system. Then, the pulsar terms can largely widen 
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Figure 1. Schematic pictures of primordial black hole (PBH) 
searches using a pulsar timing array (PTA). The impulse signal 
of a PBH is characterized by its mass MppH! the relative veloc- 
ity V to the target, the impact parameter b, and the projection 
angle 9 between the pulsar-Earth line and the closest approach. 
Top panel shows the case with a PBH passing nearby the Earth. 
The acceleration by the PBH is imprinted in all the timing data 
available. By taking the correlation of the timing data, the signal 
can be effectively amplified by a factor of V A/pSR- Bottom panel 
shows the case with a PBH passing near a pulsar. The PBH mod- 
ulates the timing data of the specific pulsar alone. The rate of such 
encounter is proportional to the number of pulsars Npsn- 

our survey volume, and our sensitivity for the direct PBH 
search would become bett er than the previous estimation 
only with the Earth terms |Seto fc Cooravl[2007r ) . In this pa- 
per, we study this issue with special attention to the differ- 
ences between the roles of two terms for the PBH search with 
PTAs, namely, the independent and numerous pulsar terms 
and the coherent Earth terms. We find that these terms 
have both advantages and disadvantages for PBH search, 
and work in a complementary way. 



2 PROBING PRIMORDIAL BLACK HOLES 
USING PULSAR TIMING ARRAY 

In this paper, we assume that our PTAs are composed by 
totally A^psr pulsars with roughly the same level of tim- 
ing noises. We further assume that the noises are white 
and have no correlation between different pulsars (see 



IShannon fc Cordesl (|20ld) for impacts of red noises). Note 
that the total numbers of the Earth terms and the pulsar 
terms are both Apsr. 

In Fig.l, we provide schematic pictures of PBH search 
with a PTA. The top panel shows the case in which a PBH 
passes nearby the solar system and excites the Earth terms 
coherently for the TOA data of all the pulsars. In the bottom 
panel a PBH fiybys a pulsar and modulates only the specific 
pulsar term. 

As mentioned earlier, we can amplify the sensitivity of 
the Earth terms by a factor of ~ 1/ViVpsR using their co- 
herent structure. In contrast, the fiybys of PBHs around the 
whole pulsars occur ~ Apsr times more frequently than 
those around the Earth alone. 

We hereafter assume that, when a PTA is available, the 
coherent Earth term can be removed from TOA data of each 
pulsar, and its pulsar term can be analyzed separately. 



2.1 Signal to noise ratio 

Here we estimate the signal-to-noise ratio (S/N) of a flyby 
event in the TOA data of a PTA. The target mass can 
be the Earth or one of the pulsars. In the Fourier space, 
the main contribution of the impulse acceleration is the 
mode with the frequency / = 1/T, where T is the time 
scale for the PBH passing around the target. In the same 
manner, the amplitude s/ of the mode can be estimated as 
Sf ~ (aT 2 /c) x | cos0|. Here a is the peak magnitude of ac- 
celeration and 9 is the angle between the pulsar-Earth line 
and the closest approach (see Fig.l). With the impact pa- 
rameter b and the relative velocity V, the two quantities T 
and a are given as 

T 4- 1 ^<74^)(35^m7i)" 1 ' « 

and a « GAfpBH/6 2 , where Mpbh is the mass of the PBH. 
Then the Fourier amplitude of the timing residual is ex- 
pressed as sf ~ (GMpbu/cV 2 ) x |cos#|, or 

1fl /Mpbh\ / V_ y 2 f\cos8\\ 

Next we evaluate the noise associated with the timing 
analysis. Using the sampling rate v of TOAs and the rms 
noise a of each TOA, the Fourier mode of the timing noise at 
the frequency / = 1/T is given as ~ a/VTu. This result is 
for TOA data of a single pulsar. When we deal with a PTA, 
it is important to distinguish whether we examine the Earth 
terms or a pulsar term. For the former, we have the statis- 
tical reduction factor 1/V5Vpsr from the coherence of the 
signals. Therefore, the effective noise level for the impulse 
search can be expressed as 

nj „ 6 . 2ns ( JL\ - 1/2 / u yV* B , 

1 VIOOns/ \Wyr J VO.SOwk^ 1 / 

(3) 

with E — or 1 for a pulsar term and the Earth terms, 
respectively. From Eqs.Q and © the signal-to- noise ratio 
of the flyby detection is now given as S/N = s//n/. 

The signal of PBH is observationally characterized by 
the duration T and the amplitude s/. On the other hand, 
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there are four physical parameters, b, Mpbh, V, and #0 
Thus, in general, we can obtain only two constraints be- 
tween these four parameters, even if the signal is detected 
with a high signal-to-noise ratio. However, when constrain- 
ing the parameters of PBHs as a dominant component of 
dark matter, we can set fiducial values for V and 6 based 
on the following considerations. The rms velocity for halo 
dark matter relative to th e solar system is dynamical ly es- 
timated to be ~ 350km/s ()Carr fc Sakellaria dou 1999J). Be- 
sides, the typical peculiar veloc ity of the observed MSPs 
are relatively small < lOOkm/s (|Hobbs et alj 120051 ) . Thus, 
it is reasonable to fix V = 350km/s for discussing both the 
Earth and pulsar terms induced by PBHs. Furthermore, if 
the scatterings between PBHs and the targets occur isotropi- 
cally, the ensemble average of the projection angle 9 becomes 
,/(cos 2 6») = l/v/3 ~ 0.58. Hereafter we fix | cos6»| = 0.58 as 
the fiducial value. Now the observational parameters (T and 
Sf) and the physical parameters (6 and Mpbh) have one-to- 
one correspondence in our order-of-magnitude estimation. 







2.2 Event rate 

The density of the lo cal dark matter is e stima ted to be 
Pdm = 0.011M a pc~ 3 (Oiling fc Merrifieldl (|200ll ). see also 
iMoni Bidin et aD (|2012 ) for a recent claim) . We put r\ as the 
mass fraction of PBHs among the Galactic dark matter, and 
assume that the PBHs have an identical mass parameterized 
by Mpbh- Then we can estimate the event rate of the close 
encounter of a PBH around the target masses within a im- 
pact parameter b as R « Tvb 2 V x (/Opbh/A/pbh) x Npsb_ 1 ~ e , 
or 



R 



0.032yr _1 
( b 



Pdm 



V740AU7 V350km/s 



O.OllM pc- 3 
V 



10 25 E 



iVpSR 



(4) 



We take E = or 1 depending on whether a PBH flybys 
one of the pulsars (E = 0) or the Earth (E = 1). The factor 
Npsn 1 ~ E reflects the fact that, only in the case of using the 
independent pulsar terms, the event rate is proportional to 
the number of pulsars. 



2.3 Detectable parameter regions 

Now we discuss the observational prospects of PBH search. 
To this end, we examine the detectable PBHs in the two 
dimensional space (Mpbh, 6), assuming two future PTAs 
whose basic parameters are summarized in Table. 1. PTAa 
has a relatively conservative se t of parameters, and could 
be re alized in the near future (|Ellis. Jenet. fc McLaughlinl 
12012} ). The goal of PTAb is more challenging and might 



1 For a PBH passing around the Earth, we can, in principle, 
estimate the direction of its closest approach, using the dipolc 
pattern of the Earth terms in a PTA. 

2 We should note that some of the observed MS Ps have com- 
parab le or larger peculiar velocities (e.g. B1957+20 iHobbs et al.1 
1120051 )). Also the scattering between PBHs and the targets may 
occur in a non-isotropic way, in which -J (cos 2 9) takes a different 
value. We discuss these cases in the final section. 



Table 1 . Parameter sets of the two pulsar timing arrays assumed 
for future prospects of PBH searches 





PTAa 


PTAb 


timing noise in each TOA (<r) 


100ns 


10ns 


frequency of TOA sampling (u) 


O.Swk" 1 


l.Owk" 1 


number of pulsars (ATpgjj.) 


100 


1000 


observation time (T a ^ a ) 


lOyr 


20yr 



be rea lized with future facilities like SKA (e.g. ISmits et al.l 
(^Ollf )). The latter will also enable us to estimate the par- 
allax distances to the stable pulsars with typical errors of 
< 20%. Then we can obtain broad scientific results^ includ 



inga map of the interstellar electron density (|Smits et al.l 

For detecting a PBH, we request that the follow- 
ing three conditions are simultaneously satisfied in the 
(Mpbh, 6)-plane; 

(i) S/N is larger than a threshold value (e.g. S/N > 3 for 
99% confidence level). 

(ii) At least one event occurs during the observation time 
(RxT ohs > 1). 

(iii) The duration of signal is shorter than the observation 
time (T<T obs ). 

From Eqs.Q and (|3}, the condition (i) corresponds to 
the region in the (Mpbh, fe)-plane as 



M P 



10 25 g J V740AU 



1/2 



> 0.56 



( V 



100ns 

5/2 



V 350km 



0.5wk _1 

| cos 9. 
0.58 

E/2 



-1/2 



(5) 



By reducing the rms noise a, the overall sensitivity is im- 
proved and smaller PBHs are within reach. With the Earth 
terms, the detectable PBH mass becomes 1/ V-^psr time 
smaller than the limit with a pulsar term. 

From Eq.Q, the condition (ii) can be expressed as 



M P 



10 25 g J V740AU 



< 0.32 



Pdm 



O.OllMgpc- 3 

/ Jobs \ 



x Apsr 



350km/s J VlOyr / 

l-E 



(6) 



When fixing the PBH fraction j), the event rate scales as 
R oc AfpBH -1 , and PBHs encounter the individual target 
masses less frequently for larger Mpbh. On the other hand, 
the event rate is also proportional to the number of the 
target masses. Thus, rare encounters with massive PBHs 
can be probed with the pulsar terms. 

From Eq.{TJ, the condition (iii) for the signal duration 
can be expressed as 



b<bn 



740AU 



Tob 



V 



(7) 



lOyr J V350km/s / 

When we fix the relative velocity, the signal duration T 
has the one-to-one correspondence to the impact parame- 
ter b(= VT). The observation time limits the duration of 
the detectable signal, and so does the impact parameter. 
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Figure 2. Detectable parameter regions of a primordial black 
hole (PBH) with mass Mpbh passing nearby the Earth (top 
panel) or a pulsar (bottom panel) with impact parameter b. 
The dark-colored regions correspond to the detectable events by 
PTAa, and the light-colored ones for PTAfo (see Table. 1). Here we 
fix the relative velocity of PBHs at V = 350km/s, and the projec- 
tion factor at | cos 8\ = 0.58. We assume that the dark matter den- 
sity is p DM = O.OllM pc _3 with the PBH fraction r) = 1. The 
lines with the slopes b oc Mpbh -2 , b oc Mpbh 1 ^ 2 , and b = const 
correspond to the conditions (i), (ii) and (ii) in the main text, 
respectively. 



Fig. 2 shows detectable parameter regions of a PBH in 
the (Mpbh, 6)-plane for r\ = 1. The top panel is for a PBH 
passing near the Earth, and the bottom one is for those near 
the pulsars. In each panel, the dark-colored regions represent 
the detectable events with PTAa, and the light-colored ones 
are for PTA6. 

In Fig.2, the lines with b oc M PBH ~ 2 , b oc A/ PB h 1/2 , 
and b — const correspond to the conditions (i), (ii) and (iii), 
respectively (see Eqs.Q, (0 and (Q). As one can see easily 
from Fig.2, the Earth terms and the pulsar terms are highly 
complemental. The Earth terms cover a lighter mass range 
~ 10 22 ~ 24 g, compared with the range ~ 10 24_28 g probed by 
the pulsar terms. 

Fig. 3 shows possible constraints on the PBH fraction 
77 with the two PTA models. Again, the top (bottom) 
panel shows the case using the Earth (pulsar) terms. The 
detectable region is bounded by two lines. The condition 
S/N> 3 sets the minimum mass of detectable PBH with the 



Figure 3. Possible constraints on the primordial black hole 
(PBH) abundance (77 = Ppbh/pdm) using the two PTAs (the 
dark-colored regions for PTAa, and the light ones for PTAfe). 
As in Fig.2, the top (bottom) panel shows the case with the 
Earth (pulsar) terms. We put V = 350km/s, |cos0| = 0.58, and 
PDM = O.OIIMqpc -3 . The lines with the slopes Mpbh = const, 
and rj oc Mpbh correspond to the conditions (i) and (ii) respec- 
tively, with the impact parameter at fe max . 



line Mpbh = const (see Eq.([5}), and event rate determines 
the maximum one with 77 oc Mpbh (Eq.©). These bound- 
aries can be understood in the following manner. For a given 
observational period T b s , the impact parameter takes the 
upper limit 6 ma x = VT h B given in Eq.©. Next, when the 
PBH mass Mpbh is fixed, the signal-to-noise ratio and the 
sensitivity to the ratio 77 both take their optimal values at 
b = 6 m ax, as shown in Eqs.([5]) and JBJ. 

If no event is detected by the PTAa observation, PBH 
with 10 24 g to 10 26 g is excluded as a dominant component 
of dark matter (77 < 0.1). Moreover, in the case of the PTAb 
observation, PBH within 10 22_28 g could be excluded (77 < 
0.01). 



3 DISCUSSIONS 

Our order-of-magnitude estimation has shown that we might 
probe PBHs or constrain their abundance with future PTAs. 
Here let us discuss issues related to the PBH search. 



Exploration for PBHs using PTAs 5 



3.1 Data analysis and competing noises 

So far we have simply evaluated the signal-to-noise ratios of 
fly-by events, assuming white noise spectra for timing noises 
of MSPs. But we should pay much attention to whether PBH 
signals can be clearly identified in their TOA data. 

The actual timing analysis typically proc e eds a s 
follows (see e.g. iHobbs. Edwards, fc Manchester! (|2006h ; 
lEdwards. Hobbs. fc Manchester! (|200^ 7 for the detail). First, 
the measured TOAs are converted to the pulse emission 
times in the presumed reference frame of each MSP. In this 
step, various effects are estimated, such as the propagation 
delays, the ephemerides, and the coordinate transformation 
from the Earth to the pulsar. The derived time of emission 
is fitted by a timing model including a pulse frequency and 
its time derivative. The residual is divided into a white noise 
and un-modeled systematic noises which could include the 
signal we are seeking. A key in our approach to enlarge the 
detectable range of PBH mass is to consult both the coher- 
ent Earth terms and the individual pulsar terms. These two 
are assumed to be separated. 

Here, we should be aware of the risk to mistake PBH sig- 
nals for other effects while processing the TOA data, which 
could also lead to misestimation of the parameters of each 
pulsar. To prevent this, more detailed modeling of the PBH 
signal is required. 

With a PTA, we have a large number of independent 
pulsar terms. The impulse accelerations by fly-by PBHs ap- 
pear in the pulsar terms as sparse and isolated signals with 
finite durations. Therefore, a long-term observational cam- 
paign would be quite helpful to select sufficiently stable 
MSPs and examine the characteristic time profiles of the 
PBH signals. 

As mentioned above, we have assumed that the ordi- 
nary timing noise in each TOA is white, own oc T° (T: 
the observational time span). As for non-recycled pulsars, 
there have been confirmed secular red noises, <tr,n oc j i2±0 - 2 
( Sha nnon fc Corded I2OI0I ) . Although most of the observed 
MSPs only have upper limits on the amplitude of the red 
noise, it might be identified by future PTAs to be obstacle 
in the PBH search. 



pattern (I — 1) whose direction is determined by the accel- 
eration vector. On the other hand, the GW background will 
have multiple modes starting from the quadrup ole (1 = 2). 

One may suspect that floating planets (|Sumi et al.l 

l201l! ) could also become an effective noise because of the 
similar masses to PBHs we are interested in. There are 
expected to be a few times more floating planets than 
stars. The anticipated encounter rate with a pulsar is ~ 
10 _(7 " 8) yr _1 x (&/740AU) 2 , which is much smaller than PBH 
constituting a dominant fraction of dark matter (see Eq.Q). 
Thus, we conclude that floating planets would not be prob- 
lematic for our approach. 

3.3 Possible refinements 

Finally let us discuss possible improvements for our estima- 
tion of the detectable region of the PBH parameters. 

In Fig. 2 and Fig. 3, we have fixed V = 350km/s, 
| cos(9| = 0.58, and pdm = O.OIIMqpc -3 . These values cor- 
respond to the observed rms velocity for halo dark matter 
relative to the solar system, the mean projection factor for 
isotropic scatterings, and the local dark matter density near 
the solar system, respectively. These treatments should be 
regarded as a zeroth order approximation, and we can refine 
the present analysis in a more realistic manner. 

For example, while the motions of the MSPs would be 
typically close to the Galactic rotation, those of PBHs would 
be more isotropic. Then the mean values for the projection 
factor cos 9 and the relative velocity V would depend on 
the Galactic position of each target MSP. In addition, the 
densities of the PBHs and the MSPs would also depend on 
the Galactic position. These would affect the estimation of 
the event rate and the strategy for data analysis. 

Future facilities like SKA could determined the distance 
of M SPs up to 13kpc with < 20% accuracy (|Smits et al.l 
1201 if ). Multiple detections using the pulsar terms could po- 
tentially give us more detailed informations of the PBHs, 
such as their density and velocity distributions in the whole 
Galaxy, which can not be reached only using the Earth 
terms. 



3.2 Competing signals 

Once a higher sensitivity as we have considered is realized, 
other possible signals also have to be opened up for discus- 
sion. They can be effective noises for PBH search and vice 
versa. 

The most promising target of future PTAs is the GW 
background from merging supermassive black holes. The 
typical magnitude of the estimated dimensionless strain 
is h ~ IP" 15 - 5 x (f /0.1yr~ 1 )~ 2/3 |jaffe fc Backed 120031 ; 
IWvithe fc Loebl 120031 ). This corresponds to the amplitude 
of the timing residual of s /jG w ~ 4.1ns x (// O.lyr -1 )" 5 / 3 
(for t he translation from h to s/^w, see e.g. IHobbs et al.l 
(2009)). which can be comparable to that induced by PBHs 
(see Eq.©). Given the effective sensitivities of the pulsar 
and the Earth terms, the latter would be more vulnerable 
to the existence of the GW background. But, using the an- 
gular patterns of the Earth terms, we can, i n principle, sep- 
arate the two competing signals, as follows (jSeto fc Cooravl 
l2007f ). A PBH signature in the Earth terms will have a dipole 
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